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ABSTRACT 


This  is  the  second  part  of  a  study  of  stochastic  processes 
generated  by  failures  and  repairs  of  components  in  a  series 
system.  In  particular,  we  consider  a  two  component  series 
system  where  failure  ot  component  1  shuts  off  component  2 
hut  not  vice  Versa.  In  hart  I,  selected  measures  of  system 
performance  were  calculated  for  general  failure  and  repair 
distr  ibut  Ions  as  ir.pli.  it  functions  of  system  availability. 
In  order  to  obtain  explicit  results,  wo  assume  in  this  part 
that  either  component  1  or  component  2  has  exponential  fail¬ 
ure  and  repair  distributions.  for  the  case  where  component 
1  has  a  general  failure  distribution,  we  obtain  counter¬ 
intuitive  results  concerning  the  long  run  average  of  system 
uptimes.  Intuitively,  it  is  best  to  shut  off  an  operating 
component  when  the  system  is  down.  However,  this  is  not 
correct  if  the  component  failure  distribution  is  1FR  (or 
more  generally  NBt'F)  and  we  wish  to  maximize  the  long  run 
average  of  system  uplines.  We  use  the  method  of  supplemen¬ 
tary  variables  to  obtain  explicit  results  for  important 
special  cases. 


ASYMPri'TTr  KKASURKS  OK  SVlTi'V  VFKFORVANt:!' 

:  ,v:.i  1'RNAl  iVl.  »-:-E|LvnNv:  Kli.lS,  li 

by 

Richard  K.  Barlow  and  Fsther  Si.!  Kudos 

1 .  INTRODUCTION  AND  SUMMARY 

In  this  paper  we  study  stochastic  processes  generated  by  failures 
and  repairs  of  components  In  a  series  system.  A  series  system  of  k 
components  operates  If  and  only  If  each  of  the  k  components  operates. 
However,  depending  on  the  shut-off  rule,  some  components  may  continue 
to  operate  with  the  system  down.  For  example,  failure  of  the  power 
supply  may  shut  down  a  computer  but  not  vice  versa.  Only  failed  com¬ 
ponents  arc  repaired  or  replaced,  and  repair  or  replacement  takes  a  random 
time.  Repaired  components  are  assumed  to  function  like  new  components. 
Furthermore,  components  are  separately  maintained.  Failure  and  repair 
times  are  statistically  independent. 

We  are  interested  In  the  asyraptot ic  (a.  time  becomes  infinite)  values 
of  selected  measures  of  system  per f ormance .  Some  o(  these  quantities 
are : 


(i)  The  limiting  system  availability;  that  is,  the  limiting 
probability  that  the  system  is  functioning; 

(11)  The  limiting  system  failure  rate; 

(ill)  The  limiting  average  of  system  uptimes  and  downtimes; 
(tv)  The  limiting  average  number  of  system  failures  due  to  a 


spec i f led  component . 


2 


Those  measures  have  various  uses.  For  example,  (Iv)  might  be  used 
C u  evaluate  the  relative  importance  oi  ditJerent  components  In  the  system. 

In  Part  1  |1],  we  computed  these  measures  implicitly,  for  systems  with  two 
components,  as  functions  of  limiting  system  availability.  See  [1)  tor 

notation  and  related  work.  In  this  part  we  assume  that  either  component  1 
or  component  2  has  exponential  failure  and  repair  distributions  and  give 
explicit  formulas  for  the  above  measures. 

The  basic  results  obtained  are  that: 

•  Asymptotic  measures  of  system  performance  defined  above  depend 
on  the  Laplace  transforms  of  probability  distributions  at 
specified  values  of  the  argument  as  well  as  means; 

•  Sharp  hounds  on  asymptotic  measures  of  system  performance  can 

bo  Riven  in  terms  of  means  assuming  certain  distributions  belong 
to  the  class  of  "new  better  (worse)  than  used  in  expectation" 
and  "new  better  (worse)  than  used."  [Marshall  and  Proschan 
(1972).  ] 

•  If  component  l  is  NBUF  and  we  want  to  maximize  the  long  run 
average  of  system  uptimes,  then  we  o)uuld  not  shut  off  1  when 
2  falls.  If  component  1  Is  NWl'E,  then  we  should  shut  it  off 
when  2  falls. 

Classes  of  life  d f s t r i but  ions  described  in  Barlow  and  Troschan  (1975) 
are  the  key  to  obtaining  bounds  on  asymptotic  measures  of  system  performance. 


A  life  distribution  F  with  mean  ti  is  said  to  be 
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(i)  Increasing  ■•ere  ..•> ing)  failure  rate  or  Irk  (DFK)  if 

log  F(x)  is  concave  (convex)  for  x  >  0  ; 

(ii)  Increas  in;.;  (decreasing)  failure  rate  ivorage  >r  !  I'RA  (DFRA) 
if  -  l log  F(x)l/x  is  Increasing  (decreasing)  in  x  >  0  ; 
(ill)  New  better  (worse)  than  used  or  NBU  (NWU )  if 

F(x  +  y)  1  (>)  F (x) F(y ) 

for  all  x  ,  y  >  0  ; 

(iv)  New  better  (worse)  than  used  in  expectation  or  NBUK  (NWUF.) 
if 


for  all  l  >  0  such  that  Fft)  >  0  . 

(v)  Discounted  life  is  less  (greater)  than  the  exponential  or 

MX  (M.C.)  if 


f  e"StdF(t)  <  (>)  ~ -  - 

J  1  ♦  su 

0 

for  all  s  >  0  . 

The  chains  of  inpl lent  tons 

IFR  «♦  1  FRA  ■—  NBU  -*  NBUF.  -*  DLL 
DFR  —  DFRA  NWU  -*■  NWUF  -*•  DLC 

are  easy  to  establish.  For  completeness  we  prove  the  last  implication 


in  each  chain. 
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l.etnma  1.1: 

If  F  is  NBUE  (NUUE)  with  mean  y  ,  then 

o> 

/  i  <i>  i-fj: 

0 

for  all  s>0. 

Proof : 

For  F  NBUF, 


yF(t) 


F( t  +  x)dx 


/ 


F(x)dx  . 


C 


Hence , 


u 


f e  StF(t)dt 
0 


OB 

/ 


F(x)dxdt  . 


t 


Interchanging  order  of  Integration  on  the  right  and  reducing  yields 


<a» 

/  0 

0 

"RtF(t)dt  >  -j— —  . 

-  H  us 

OD 

But  f  e  S?F(t)dl  ■ 

0 

j^l  1  dF(t  )j  /s 

so 

m 

that  f  e  StdF(t)  < 
0 

1 

l  +  us 

All  inequalities  are  reversed  if  F  is  NWUE.H 


with 


We  assume  that  component  i  has  failure  distribution  F. 
mean  -  1  /  a  ^  and  repair  distiibution  t;^  with  mean  -  l/o^  . 

I  “  1,2  .  Let  1'j,L'2 . U  be  successive  system  uptimes  under  our 

policy  tluit  1  shuts  off  1  but  not  vice  versa;  i.e.. 


o  o 


1  2 


Let 


1  ira 
n 


n 


+  U 

n 


when  this  limit  exists  almost  surely.  In  Tart  1  |1),  we  showed  that 
when  F  is  exponential 

1 

Li  ■■  -  . 

•*  "l  +  *2 

This  is  also  true  tor  general  distributions  and  the  shut-off  policy 

o  o 

1  2 

i.e.,  neither  component  shuts  off  the  other,  as  well  as  the  shut-off  policy 

o  -  —  o 

1  2 

I.e.,  failure  of  either  component  shuts  off  the  other. 

A  major  result  of  this  paper  is  that  if  Fj  is  NBUF  (NWUE)  and 
F,  ,  r.  are  exponential,  then 
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Actually,  flic  Inequalities  1 J  for  1  ^  DLL  (DI„G)  .  IIimkc t  in  terms 

of  limiting  average  of  system  uptimes,  it  is  better  not  to  shut  off  1 
when  2  fails  and  Fj  Is  NBUE. 

In  Part  I  (1),  we  showed  that  in  general 


* 


/> 

\  1  Vo 


Theorem  2.1  provides  explicit  expressions  for  and  when  F^ 

and  Gy  are  exponential.  figure  1.1  is  a  graph  of  \i  as  a  function 
of  a  where  F ^  has  the  gamma  density 


f  r  (x) 


a  a-1  -ax 
ax  e 


r<  >) 


Note  that  F  has  mean  1.  The  graph  is  for  X  -  X  •  p  -  ~  i 

1  12  12 

And  cj  •  r?  •  ^2  exponential.  F(  in  this  case  is  IFR  (and  hence, 
NBl'E)  for  a  >  1  .  For  0  <  a  <  1  ,  F  is  DFR  and  hence,  NWUE. 

Figure  1.2  is  a  graph  of  limiting  availability  A  ^  as  a  function  of  a 
umder  the  sane  assumptions. 

Of  course,  other  considerations  such  as  the  cost  of  repairing 

component  1  may  make  our  solution  less  attractive.  In  general,  the 

limiting  availability  seems  always  to  be  greatest  for  the  policy  o  -  —  o 

1  2 


% 


KX!’  'NKN  !'  I  A’. 


2.  F, 


In  this  section,  we  assume  that  both  failure  and  repair  distributions 

iponent  2  ar«  mantial,  and  ,  .  cel 

with  failure  rates  l^(t)  ,  6^(t)  ,  respectively.  We  obtain  explicit 
expressions  for  the  stationary  state  probabilities  In  general,  and 
for  the  limiting  availability  In  particular. 

We  obtain  bounds  on  the  availability  A(FpC^)  for  our  model.  We 
also  compare  the  other  stationary  probabilities  as  well  as  u  and  v 
with  the  corresponding  quantities  for  the  all-exponential  case. 

Finally,  we  compare  the  quantities  A  ,  p  ,  v  with  the  corresponding 
quantities  in  Models  A  and  h  of  Fart  1.  We  use  the  method  of  supplemen¬ 
tary  variables  described  In  Cox  (19SS)  to  obtain  stationary  probabilities. 

Following  the  notation  in  I'art  1,  * ^  is  the  limiting  probability 
that  component  1  (i  “  1,?)  alone  Is  under  repair,  is  the 

probability  no  component  is  under  repair  and  i  ^  is  the  probability 
both  components  are  under  repair.  See  Fart  1  for  additional  notation. 


r  “St  - 

We  need  the  notation  F.(s)  "  J  e  l(t)dt  - 

1  0  * 


1  -  f,(s) 


where 


f  .  ( s )  »  f  e  h  d  F  ( t  )  . 


Theorem  2 . 1 


F  ,  ,  G,  are  exponential  and  Fj(t)  «  e 


-/  1 » (u)du 

0 


-/  0, (u)du 
-  0  1 

Gj(t)  ■  e  ,  then  the  stationary  state  probabll ities  are  given  by 


i  ind  the  onvenrion 


p(t,u)  *  0  i  i>r  u  t)  and  /  p0(t,l^)dt^  «• 

S  p  (t,S  )d*  ♦  /  p,(t,t  )dt  ♦  5  p.ft.s  )ds  -  1  . 

0  1  0  0 

Since  we  are  only  Interested  In  the  stationary  probabilities,  assuming 
they  exist,  we  note  that 

5p  (t,u) 

lim  — ^-r- -  "  0  . 

Q  t 
t 


Denote 


p  (u) 
rn 


def 

H 


i  la  p 


(t  , u ) 


and 


n 


p  (u)du 
n 


the  stationary  probability  ot  being  in  state  n  ,  n  “  0,1,2, 1  . 

letting  l  *  •  ,  we  get  the  differential  equations  for  the  stationary 
Joint  probabilities 

1P^(u)  -  -(Xjlu)  +  >,)p()(u)  +  e2p2(u) 
pj(u)  -  -Aj (w)Pj (u)  +  O^p j  (u> 

Kt.H)  \ 

j  p',  (n)  -  XjP^ln)  -  (Xj(u)  +  02>p,(u) 


P>) 


-(e  (u)  +  e_)p-(vi) 


1  3 


ilh  t  he  initial  lit  i  .»ns 


(2.12) 


/ 


I  ,  : 


Pj(0)  -  J  1 i (u)pQ(u)du 


8) 

>2(0)  ’J  0! 


(u)pj(u)du 


P  JO)  "  J  J  j  (u)p2 
0 


(u)  dll 


and  the  normalizing  equation 


”0  +  T1  +  *2  +  '3  -  1 


Solving  the  last  equation  In  (2.11)  yields 


P3(t)  -  BjC^Oe 


-M 


(2. 1 3) 


OB 

>  3 (0)  -  B3  -  J  (u)p2(u)dvi 


from  the  last  equation  in  (2.12). 

Assuming  for  the  moment  Oj(u)  0^  ,  i.e.,  constant  repair  rate, 
the  characteristic  polynomial  of  the  second  and  fourth  equations  in  (2 

pj  -  -V,  +  v3 

P3  ’  -(61  +  62)p3 
Is 


(  2 . 1 A  ) 


.11) 


( x  +  0  )  (x  0  l  f  1  , ) 


with  roots  -0  ,  -(0,  +  09)  .  so  that 


(2.15) 


-*  ,t  -(0  +0.,)t 

l>  1  ( t- )  “  A  j  e  +Bje 


-  Cj(t)  Al  +  BjC 


This  is  easily  seen  to  be  the  solution  in  the  case  6^(u)  t  6^ 

(with  the  help  of  (2.13)),  and  substituting  back  yields  conditions  on 

A1  *  B1  *  B3  = 


(2.16) 


P3  "  "B1  ■ 


The  first  and  third  equations  in  (2.11),  lor  A^(u)  ,  are 


(2.17) 


(  P0  ’  "(X1  +  VP0  +  °2p2 

(  p2  “  Vo  '  (X1  +  °2)p2 


with  the  characteristic  polynomial 


( x  ♦  A j  +  A7)(x  +  Aj  +  0^)  -  AjOj  "  (x  +  Aj)(x  +  ),j  +  Xj  +  0^) 


Hence,  the  general  solution  of  (2.17)  is 


(2.18) 


-At  -(A  +X,+e.)t 

p  (t)  -  A  e  +  Be 

n  n  n 

T  ~(X2+B2^tl 

*  F ,  ( t )  I  A  +  B  e  I ,  n  •  0 , 1 


Substituting  Kk  \  in  (2.11), 


'  10  .  e  :  l  si-eii  to  he  a  1  so  the  solution  for 


X^(u)  t  Aj  ,  with  the  conditions  on  the  coefficients 


(2.19) 


B0  +  B2  ‘  ° 


'Vo  +  °2A2  “  ° 


Substituting  (2.13),  (2.15),  (2. IS)  in  the  initial  conditions,  Equation 
(2.12)  yields  the  following  relations  among  the  constants 


\)  4  *0  “  A1  +  BlW 


(2.20) 


A1  4  B1  “  A0  +  Vl(X2  4  V 

a2  +  b2  -  b,r*(o2) 


B3  “  A2  4  B2fl°2  4  V  ’ 


Integrating  (2.13),  (2.15),  (2.1R)  gives  four  equations  for  the  n  's  , 

n 

and  together  with  the  normalizing  condition,  we  have 


*0  ■  i;*  V'>2  *  V 


'L  *  *  »,V°2> 


(2.21) 


■»  *  * V 


•3  '  WV 


Equations  (2.19),  (2.20),  (2.21)  form  a  system  of  eleven  linearly 


independent  equations  in  the  eleven  unknowns 
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VV'i’ VAo,Bo,Ai,urAr82,B3  ‘ 


The  solution  of  (2.19),  (2.20),  (2.21)  for  the  unknown  coefficients  yields 


(2.22)  < 


1  -  »;<>2  +  »;>«>2>  . 

A0  -*  Bo 

hW 


_*  ®Q 

WV 


Ai  -  o  +  p2)  — r* 


X2Gl(02) 


A2  ■  — 


wv 


K  (X  ♦  0  ) 

Hi"(l  +  P2>'T>7"Bo 


B2  "  "B0 


F  (X  +  0  ) 

B  -  (1  +  p  )  -I-; - B 

Cl(92> 


and  substituting  (2.22)  in  (2.21)  yields  the  expressions  for  v  ’s  in  the 

n 


theorem,  as  well  as 


(2.21) 


„  _ 

0  S  ♦  0!  X2  ♦  02 


*2*2  VV 


To  get  tj,  ,  itj,,  ,  split  the  equations  dealing  with  Pj(u)  in  (2.11), 


(2.12)  into 
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pj  (u)  -0  (u)|>  ,  (u) 


Pllt(u)  -  -6j (u)p1 , , (u)  +  e2p3(u) 


(2.24) 


JX  l(u)pQ( 


P[ I (0)  -I  Aj (U)p  (u)du 


so  that 


pp  ,  (0)  -  0 


Pj.  (O  -  Aj.GjCt) 


Pj,(0)  -  Aj,  " (u)pQ(u)du 


“  A0  +  Vl(V2  +  °2> 


Ar  Ao  .  Bo  *, 


r  "  8'1  "  Oj  +  p~  ‘l(v2  +  e2J  ’ 


which,  using  (2.22),  (2.23),  can  be  shown  to  be  as  stated  in  the  theorem. 
Also,  similarly  to  (2.13),  we  get 


-  G,<0  A,..  ♦  »,..<• 


’[*1" 


0  -  Aj , ,  +  Bj , , 


Substituting  back  in  the  second  equation  of  (2.24)  and  using  (2.13)  yields 


V’  "  'B3 


so  that 


A1 ' •  _* 

-  -07-  *  B1--C.t02> 


“i 

-O’  11  -  WV 


which,  using  (2.22),  (2.23)  can  bt*  shown  to  be  as  stated  in  the  theorem. 


Remark : 

From  the  preceding  theorem,  we  can  find  explicit  expressions  for 
the  p^(u)  •  n  “  0,1, 2, 3  (and  n  -  1  *  ,  1  *  * )  .  That  is,  the  stationary 
probabilities  of  being  in  state  n  ,  and  the  age  of  component  1  being 
u  ,  for  n  "  0,2  ,  and  the  stationary  probabilities  of  being  in  state 
n  ,  and  the  repair  time  of  component  1  is  u  ,  for  n  “  1* ,1,,,1,3  . 

This  can  be  achieved  easily  by  substituting  the  values  of  the  constants 
from  (2.22),  (2.23)  in  the  expressions  for  p^u)  (2.13),  (2.15),  (2.18). 

In  Part  1,  wo  showed  in  Theorem  2.2.10  that  for  arbitrary 
distrlbut  ions 


(2.25) 


1 _ 

°lV 

X  +  X, 

1  Vo  2 


Also,  it  is  known  that 


u 


for  arbitrary  distributions.  Using  Theorem  2.1,  we  can  prove 


Tiic  or  or.  .  J  : 


If  Kj  is  NBUE  (NWUF.) ,  then,  for  arbitrary  Cj 


(2.26) 


Proof : 


U  >  (<) 


\  +  X2 


err 

From  (2.25),  we  see  that  we  need  only  show  that  r— - —  1  (l)  1 

Vo 

From  Theorem  2.1, 


'lV  d  +  X2C1<02)  1 1  "  (x2  +  02)?1(X2  +  02)1 
To  d  +  VVV  VVX2  +  V 


If  and  onlv  if 


Vh  ♦  v  '  + .  . .  rt 


The  last  inequality  is 


true  for  F^  NBUK  by  lemma  1.1.1 


Remark : 


From  Theorem  2.1,  it  is  easy  to  verify  that  limit  in*  system  availability 
•  ^ ( F j , C j )  can  be  expressed  as 


AU^.Gj)  - 


(2.27) 


(1+Oj)  Yl+o, 


vw* 


e,  ..  /  e,  .. 

®7  -*  /  X2  \-*  -* 

» *  r  k:  C1  <,2>-(v*i  o.tt.P'i  <  WC1  <V 


Lemma  2  .  3 : 


for  .inv  Cj  with  mean  Vj  and  Is  the  degenerate 

distribution  with  mean  Uj  .  The  Inequality  is  sharp. 

Proof : 

(?.?8)  Is  sharp  since  we  have  equality  when  fj  Is  exponential. 
The  Inequality  follows  from  l.emm.i  2.3. 


i  .  '  •  .  •  : 


1  K !  ‘  \  in  .1  limiting  sense?. 


A l  .;>>  , 


Fj  I  FRA  implies 


is 


r  ,  1 

I  -  .  .  -sx ,  1  -  e 

J  — ; — • 

0 


Theorem  2.5: 


-1 


(i)  If  is  NBUK  (NWUE)  with  mean  “  °j  • 


then 


A ( F j .Gj)  >  (<)  (1  + Oj )  \1  + P, 


°2  /  G2  \-* 

V  •  *  ( 'i  •|iyijri<>i*li> 


for  any  F^  with  mean  l>  ^  .  The  Inequalities  are  sharp, 
(ii)  If  C  is  1  FRA  with  mean  ,  then 

AIF|.C.,>  i  a(f,."vJ 

for  anv  F  with  mean  u.  and  I)  is  the  degenerate 

1  Vj 

distribution  with  mean  Vj  .  The  Inequalities  are  sharp. 
Theorem  2_.  s : 

If  both  F  and  Gj  are  NRl!F,  then 


(2.  11) 


Tn(«,^r) 


where  »  (exp)  means  that  probabilities  are  computed  in  the  all 
exponential  ease. 

The  inequalities  are  reversed  ii  both  K-,  and  C,  are  NWUE 
The  proof  is  straightforward  and  is  omitted. 


Also, 


(3.6) 


v*  ■  piW\  f  vno 


when*  C 


(V  +  3  )  ■  f  e 

,v  1  r  *i 


“  “(Xl+61)t 


C.,(t)dr  . 


Tl>e  proof  will  he  omitted  since  li  is  similar  to  that  of  Theorem  2.1 


Theorem  ).2: 

(i)  If  G,  Is  NBL’F.  vitli  r"jn  v  7  ,  then 


(1.7) 


(l  *  p.Xl  0-)  •’ 


A(F-.G-)  <  A  (exp) 


If  G^  Is  1  FR.\  with  nean  v7  ,  then 


A(F  ,.(•,)  *  a/g,  -  D  \ 

21  l  2  v;/ 


111)  If  G,  is  NVt'F  with  mean  v,  then 


(3.8) 


i  A  °  • 1 1 2 '  ‘  A(**P) 


S’ roof  : 

TS\e  left-hand  bounds  In  (1.7)  and  (1.8)  follow  from  Fart  I, 
Corollaries  2.2.13  and  2.4.2  respectively.  They  are  always  valid  but 
not  neccssirily  sharp. 

All  other  Inequalities  are  sharp  and  the  proofs  are  similar  to 


Theorem  2.  •  .  ■ 


Theore t,  I: 


n  g2 

is  NBUK  witti  mean 

,  t  hen 

(3.‘>) 

"o  1 

"qU'XP) 

(3. 10) 

*1  i 

* j (exp) 

(3.11) 

"2  - 

(exp) 

(3.12) 

*  3  - 

*  j (exp ) 

(3.13) 

"  1  ’  - 

"j t  (exp) 

(3.14) 

•ju  i 

r  j . . (exp 

Ml  Inequalities  are  reversed  if  C ,  is  NWUF.  with  mean 


The  proof  is  straightforward  and  ir.  omitted. 
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